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ABSTRACT 


Obtaining  a  confidence  interval  for  a  parameter  A.  of  an  exponential 
distribution  is  a  frequent  occurrence  in  life  testing  problems.   Often- 
times the  test  plan  used  is  one  in  which  all  the  observations  are  censored 
at  the  same  time  point  t  .   Several  approximate  confidence  interval  pro- 
cedures are  available  in  the  statistical  literature;  however,  to  the 
knowledge  of  the  author,  the  performance  characteristics  of  the  various 
approximations  used  in  these  procedures  have  not  been  established 
analytically.   The  purpose  of  this  paper  is  to  report  the  results  of  an 
empirical  stucy  of  the  performance  of  four  of  these  procedures  with 
respect  to  the  expected  length  of  the  interval,  the  variance  of  the 
interval  length,  and  the  coverage  probability. 
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I.   INTRODUCTION 

In  life  testing  applications,  it  is  frequently  desired  to  obtain  a 
confidence  interval  for  the  parameter  A.  of  an  exponential  distribution. 
In  case  a  test  plan  is  used  for  which  all  the  observations  are  censored 
at  the  same  time  point  t  ,  several  approximate  confidence  interval  pro- 
cedures are  available  in  the  statistical  literature.   To  the  knowledge 
of  the  author,  the  goodness  of  the  various  approximations  used  in  these 
procedures,  and  hence  the  goodness  of  the  procedures  themselves,  have 
not  been  established  analytically.   The  purpose  of  this  paper  is  to 
report  the  results  of  an  empirical  study  of  the  performances  of  four  of 
these  procedures  with  respect  to  the  expected  length  of  the  interval, 
the  variance  of  the  interval  length  and  the  coverage  probability. 

The  general  setting  of  the  problem  is  as  follows:   suppose  the  random 
variables  T, ,  T„ ,  ...,  T  are  independent  and  identically  distributed  with 
common  exponential  distributions 

fi(ti}  =  e"Ui'  li  >   °;  i  =  1}  2'  ■••»  n' 

The  random  variables  that  will  actually  be  observed  are  X,,  X~  ,  ...,  X  , 

12  n 

where 


fT.       i 


f     T.   <  t 
X.    =  x  ° 


t        if     T     >  t    ;    i   =   1,    2,    ...,    n. 
o  1  o 

These   random  variables   are    sometimes   referred   to  as    "censored"   exponential 

variables,   and   t      is   called   the    "censoring   point".      Define    the   random 
o 

variables  Y, ,  Y„ ,    Y     by 

12  n     J 


>   n 


It  is  easily  seen  that  the  random  variables  Y.  ,  Y„ ,  ....  Y  are  independent 
J  1    2        n         r 

and  have  Bernoulli  distributions  with  parameter  p  =  1-e   o.   In  Section  II, 
four  confidence  interval  procedures  for  X   based  on  the  X.  and  Y.  are 
discussed . 


II.      CONFIDENCE    INTERVAL  PROCEDURES 

A.      PROCEDURE    1 

This  procedure  is  obtained  as  a  special  case  of  a  solution  to  a  more 
general  problem  that  has  been  derived  by  Halperin  [1].   The  random 
variable  Y  =  XY.  has  a  binomial  distribution  with  unknown  parameter  p. 
Suppose  y  is  the  observed  value  of  the  random  variable  Y,  then  the  100 
(l-a)$  confidence  bounds  p  and  p  may  be  obtained  by  solving  the 

U  Li 

equations 

y 


P  [Y  <  y]  -  Z  (£)  p1  (l-p)n_i  =  a/2 


i=o 
and 

P  [Y>y]  =  S  Qp1  (1-P)n]i  -a/2 
i=y 

for  p.   Due  to  the  discreteness  of  the  random  variable  Y,  these  equations 
do  not  yield  an  exact  100(l-cc)$  confidence  interval.   The  confidence 
coefficient  is  at  least  100(1-0.)$.   The  confidence  interval  for  \    is 
easily  computed  by  inverting  the  relation  p  =  1-e   o  to  get 
Xh   =  -ln(l-pL)/to  and  ^  =  -lntt-pj,)/^. 

A  useful  tool  in  solving  for  p  and  p   is  the  Incomplete  Beta 

Li  U 

function, 

F(*;a,B)  =i2±f±7li    /  ^  d-0Pdt  -1  -    S   (^f1)^  (l.x)^+H. 

^  i=o 

o 

The  upper  confidence  bound,  p  ,  is  obtained  by  solving 

F(Pu;  y,n-y-l)  =  l-a/2  ; 
and  the  lower  bound,  p  ,  is  obtained  from 

Lj 

F(pL;  y,n-y-l)  =  a/2. 


The  simplicity  of  the  computation  involved  in  this  method  is  very 
apparent;  however,  it  should  be  noted  that  the  observed  values  of  the 
random  variables  X..,  X?,  ...,  X  have  been  disregarded.   Intuitively 
it  would  appear  that  this  procedure  will  not  produce  as  good  results 
as  one  which  incorporates  all  the  information  of  the  experiment.   How- 
ever, due  to  its  simplicity,  this  procedure  merits  consideration. 


B.   PROCEDURE  2 

For  X  t  <  <  1,  Y  =  TY .    is  nearly  a  Poisson  random  variable  with 

1 


A     Y 

parameter  XIX..   Rubenstein  [2]  shows  that  X  =  ycr- 


1  +  =r 

2n. 


is  a 


nearly  unbiased  estimator  for  X,  where  the  second  factor  is  used  to 
correct  for  the  bias  of  the  otherwise  minimum  variance  estimator 
Y/XX..   Wilks  [3]  states  that  for  a  Poisson  random  variable  P  with 
parameter  p,  the  following  relationship  is  approximately  true: 
p-P  =  Zvp  where  Z  is  the  standard  normal  random  variable.   Rubenstein 
applies  this  reasoning  to  the  random  variable  Y  along  with  his  modified 
estimator  X    to  obtain  the  relationship  X  -  X   =  ZvX  .   This  equation  is 
solved  for  X   and  X   by  using  the  appropriate  standard  normal  percentage 

U  Li 


point  corresponding  to  a  100(l-ct)$  confidence  level.   The  resulting 

1/2- 


solutions  for  \      and  \      are 

L»  U 


and 


\- 


LU 


2^  +  Z2C  -  (4^Z2C  +  Z4C2) 


n 


1/2 

2t  +   Z2C  +  (4^Z2C  +  Z4C2)    H   where  C  -  (£X.)1/2 


C.   PROCEDURE  3 

This  procedure  is  a  modification  of  one  which  has  been  suggested  by 
Birnbaum  [4]  for  application  to  uncensored  life  testing  problems.   Using 
the  terminology  commonly  used  in  the  literature  of  life  testing  to 
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facilitate  the  explanation,  imagine  that  the  random  variables  X  ,  X9  , 

.  ..,  X  are  observed  sequentially.   That  is,  a  randomly  selected  item 

is  put  on  test  and  is  replaced  with  a  similar  item  at  failure  or  after 

a  period  of  time  t   has  elapsed,  whichever  occurs  first.   Then  the 
r  o 

arrival  process  of  failures  is  a  Poisson  process.   Suppose  the  experi- 
ment is  continued  until  a  total  of  k  failures  have  been  obtained.   It 
is  well  known  that  the  individual  inter-arrival  times  of  failures  are 
exponentially  distributed  and  the  time  to  the  k   failure  has  a  gamma 
distribution.   Consider  a  test  plan  in  which  the  experiment  is  stopped 
after  a  random  amount  of  time  (as  would  be  the  case  for  example,  if  n 
items  were  to  be  tested).   The  number  K  of  failures  is  then  a  random 
variable.   However,  it  would  appear  that,  given  K=k,  the  distribution 
of  the  time,  W,  ,  until  k  failures  have  arrived  can  be  approximated  by  a 

gamma  distribution.   More  precisely,  suppose  the  observed  value  of  the 

k  k 

random  variable  K  =  £  Y.  is  k,  and  let  W,  =  S  X..   Then  the  con- 

i=l  i=l 

ditional  distribution  of  W,  ,  given  that  K=k,  is  approximately  the  gamma 

distribution 

c(   N   _X k-1  -\w 

f  (w)  -  YKS  w   e   ;  w  >  o. 

It  follows  that  V  =  2\W  can  be  approximated  by  a  Chi-square  variable 

K. 

2  2 

with  2k  degrees  of  freedom.   Thus,  if  x  l-Ct/2  and  x  Oc/2  are  the  upper 

and  lower  a/2  percentage  points  of  the  Chi-square  distribution  with  2k 

2       2 

degrees  of  freedom,  then  \~~yT, >  — oTp — )  constitutes  an  approximate 

k        k 
100(l-a)$  confidence  interval  for  X. 

Birnbaum  also  suggests  an  estimator  for  X   which  is  merely  the  average 

of  the  upper  and  lower  confidence  bounds.   This  estimator  has  also  been 

computed  and  tabulated. 


D.   PROCEDURE  4 

In  Procedure  3  the  distribution  of  the  random  variable  W,  ,  the  waiting 

K 

time  until  the  k   failure,  is  approximated  by  a  gamma  distribution. 

Since  the  test  is  terminated  after  n  items  are  tested,  the  maximum  value 

that  W,  can  attain  is  nt  .   Consider  the  conditional  probability 
k  o 

P  [W <  w  and  W  <  nt  ]    P  [W  <  w] 

P  [Wk  <  w|wk  <  nto]  =  p  [u  <  „£  ]    °  =  PtH*;,tl  ■ 

k.  — •   o  k  —   o 

This  yields  a  truncated  gamma  distribution  having  density  function 

f  f(w) 

t,/\\      o  <  w  <  nt 
F(nt  )    —   —   o 

o       w  <  o,  w  >  nt 

o 

Where  f(w)  is  the  density  function  of  the  gamma  distribution  and 

nt 

r° 

F(nt  )  =  /    f(w)  dw.   Intuitively,  this  new  function  would  seem  to 

o 
approximate  the  probability  distribution  of  W  in  a  censored  test.   An 

obvious  drawback  of  this  method  is  the  difficulty  of  computation  since 
both  X  and  F(nt  )  which  depends  on  X  are  unknown.   Thus,  given  tables 
of  the  Chi-square  distribution  function  one  is  forced  to  use  an  itera- 
tive method  for  obtaining  the  desired  X  such  that 

P  [X2  <  2Xw,  ]  P  [X2  <  2Xw  ] 

■ 5 =  a/2  2 =  1_a/2 

P  [X  <  2Xnt  ]  P  [X  <  2Xnt  ] 

o  o 

to  obtain  X   and  X   respectively. 

The  estimator  for  X  described  in  Procedure  3  is  applicable  to  this 

procedure  and  is  also  tabulated. 
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III.   COMPARISON  OF  PROCEDURES 

A  Monte  Carlo  study  has  been  made  to  compare  Procedures  1,  2,  and  3 
described  above.   One  thousand  samples  of  size  n  (n  =  30,  40,  50)  from 
an  exponential  distribution  with  parameter  X    (K   =    .1,  .2,  .8,  3,  5,  10) 
have  been  generated.   For  each  sample,  confidence  intervals  for  X    have 
been  obtained  by  using  each  procedure  for  various  censoring  points  t 
and  confidence  coefficients  (l-cc)  (a  =  .05,  .01).   The  tabulated 
quantities  are  the  average  length  of  the  confidence  interval,  the  sample 
variance  of  these  lengths,  and  the  empirical  coverage  probability  of  the 
intervals  (i.e.,  the  proportion  of  intervals  which  actually  covered  X). 

An  abbreviated  Monte  Carlo  study  has  been  made  to  compare  Procedures 
3  and  4.   The  lengthy  computation  of  Procedure  4  required  that  the  number 
of  repetitions  be  reduced  to  100  and  that  fewer  combinations  of  X ,  t  ,  n, 
and  a  be  used. 
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IV.   CONCLUSIONS 

The  results  of  this  study  differentiate  each  procedure  as  to  its 
merits  and  shortcomings  with  respect  to  certain  life  testing  situations. 
Each  method  is  discussed  below  in  order  to  define  the  situations  in  which 
it  could  be  used.   The  following  discussion  includes  only  comparisons  for 
Procedures  1,  2,  and  3.   Procedures  3  and  4  are  compared  separately  and 
to  a  lesser  degree  due  to  the  differences  in  the  Monte  Carlo  studies  made 

Procedure  1  performed  as  expected;  it  generally  gave  less  accurate 
results  than  the  other  procedures  with  respect  to  all  three  quantities- 
expected  interval  length,  sample  variance  of  this  length,  and  empirical 
coverage  probability.   However,  since  the  empirical  coverage  probability 
of  this  procedure  tends  to  be  conservative  and  since  the  computation 
required  for  this  procedure  is  minimal,  Procedure  1  would  be  favored  in 
cases  whet   a  quick  but  dependable  confidence  interval  is  needed  or  when 
a  rough  estimate  is  needed  for  use  in  more  sophisticated  procedures. 
Disregarding  the  observed  value  of  the  random  variable  2X .  is  the  main 
reason  for  the  conservative  results. 

The  overall  performance  of  Procedure  2  seems  to  rank  it  first  among 
those  studied.   However,  when  t   is  near  1/A.  Procedure  3  appears  to  give 
comparable  results  with  less  computation. 

Procedure  3  performs  very  well  in  cases  where  t   is  approximately 
equal  to  l/X  ;  however,  when  t  <  <  1/A.  few  failures  tend  to  occur  and 
the  random  interspersing  of  censored  times  causes  the  sample  variance 
of  the  procedure  to  be  high.   In  general,  the  empirical  coverage  prob- 
ability for  this  procedure  is  close  to  the  chosen  confidence  coefficient. 
The  actual  computations  needed  are  comparable  to  those  for  procedure  1. 
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The  estimators  for  \   defined  under  Procedures  2  and  3  have  been 
tabulated;  the  one  given  under  Procedure  2  seems  to  be  a  nearly  unbiased 
estimator  for  \   as  stated.   The  estimator  arising  from  Procedure  3  is 
consistently  greater  than  the  true  value  of  X. 

The  comparison  of  Procedures  3  and  4  are  based  on  an  abbreviated 
Monte  Carlo  study.   This  is  due  to  the  accuracy  limitations  of  the 
computer  when  computing  small  values  of  A.  and  also  the  added  dependence 
upon  the  time  till  k   failure  in  the  iterative  step  of  the  computation. 
A  more  extensive  study  will  be  necessary  to  obtain  more  meaningful 
results.   It  appears  that  the  interval  length  obtained  from  using  Pro- 
cedure 4  is  not  as  good  as  that  of  Procedure  3;  however,  the  empirical 
coverage  probability  for  Procedure  4  seems  to  be  close  to  the  desired 
confidence  coefficient  and  the  estimator  for  \  appears  to  provide  a 
nearly  unbiased  estimator. 
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V.   COMPUTATIONAL  PROCEDURES 

The  IBM/System  360  Model  67  computer  with  the  Fortran  IV  programming 
language  has  been  used  for  all  computations. 

The  exponentially  distributed  random  sample  is  obtained  by  first 
generating  random  numbers  from  a  uniform  distribution  on  the  interval 
(0,  1)  and  then  by  inverting  these  numbers  to  get  a  random  sample  from 
an  exponential  distribution  with  parameter  A.  .   This  is  done  by  using  the 
fact  that  if  a  random  variable  Z  has  an  exponential  distribution  with 
distribution  function 


f(z)  =  p[z  <  z]  =  ; 


0  z  <  o 
l-e~Xz   o  <  z  <  1 

1  z  >  1 


the  rando  variable  X  =  1-e    has  a  uniform  distribution  on  the  interval 

(0,  1).   Thus,  if  x  is  a  random  number  from  the  interval  (0,  1),  the 

number  z  =  -  In  (l-x)A  is  a  random  sample  number  from  an  exponential 

distribution  with  parameter  \.      The  subroutine  RANDU  of  the  IBM 

Scientific  Library  is  used  to  generate  these  random  numbers. 

The  second  step  in  the  computing  procedure  requires  a  simulation  of 

censored  testing.   Therefore,  each  random  number  z  is  compared  with  the 

pre-determined  censoring  time  t  and  a  counter  is  used  to  obtain  the 

number  of  failures  k  =  £Y. .   For  Procedure  1,  only  the  value  of  k  is 

needed;  however,  Procedures  2,  3,  and  4  require  additional  information. 

If  the  value  of  z  exceeds  the  value  of  t  ,  it  is  disregarded  and  replaced 

by  t  :  these  values  of  t  and  z  which  are  less  than  t   are  summed  to 
J      o  o  o 

obtain  the  observed  value  of  the  random  variable  TX .    for  use  in  Procedure 

l 

2.   The  necessity  of  order  in  Procedures  3  and  4  require  the  random 
numbers  to  be  dimensioned  in  an  array  to  maintain  the  order  in  which 
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they  are  generated.   This  array  contains  a  random  ordering  of  z's  and 

t  *s  which  are  summed  to  the  point  of  the  k   failure  after  the  entire 
o 

sample  of  size  n  has  been  generated.   This  sum  contains  some  of  the 

truncated  times  t   but  generally  not  all  and  it  contains  all  z  <  t  . 
o  o 

The  value  of  this  sum  is  the  observed  value  of  the  random  variable  W. 

k 

of  Procedures  3  and  4. 

The  computation  of  Procedure  2  is  trivial  after  the  required  normal 
deviate  is  read  from  a  data  card.   However,  due  to  the  high  number  of 
possibilities  for  the  value  of  k  (0  <  k  <  n)  subroutines  of  the  Beta 
and  Chi-square  distribution  functions  are  used.   These  subroutines  are 
written  as  cumulative  distribution  functions  and  give  only  the  resulting 
probability,  given  the  required  input  parameter.   Since  the  input  parameter 
X    is  unknown,  it  has  been  required  to  write  an  additional  function  sub- 
program to  iterate  toward  the  desired  parameter  given  the  confidence 
coefficient.   The  restrictions  on  the  input  parameters  for  these  sub- 
routines force  the  cases  of  k  =  0  and  k  =  n  to  be  ignored.   This  is  done 
by  disregarding  the  sample  which  produces  k  and  by  accounting  for  it  in  the 
value  of  the  number  of  repetitions.   The  IBM  Scientific  Library  subroutines 
BDTR  and  CDTR  along  with  the  function  subprograms  written  by  the  author  are 
used  for  these  computations. 

For  each  of  the  thousand  samples  the  confidence  interval  length  is  com- 
puted and  is  tested  to  see  if  it  actually  covers  X.      The  final  number  of 
those  covering  X    is  divided  by  the  number  of  repetitions  yielding  the 
empirical  coverage  probability.   The  sample  variance  is  computed  by  using 
the  relation  a  =  (  2  (x-x.)  V(m-l)  where  x  is  the  average  length  of  the 
intervals,  x.  is  the  length  of  the  i    interval  and  m  is  the  number  of 
repetitions . 
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TabJe  1 


CONFIDENCE  COEFFICIENT   l-o  -  .95 


o 

PROC. 

30 

40 

50 

AVG. 

VAR. 

C.P. 

AVG. 

VAR. 

C.P. 

AVG. 

1 

VAR. 

C.P. 

3 

: 
1. 

.159 

.001 

.988 

.136 

.000 

.971 

.120 

i 
> 

.000 

.978 

2- 

.148 

.001 

.963 

.127 

.000 

.959 

.112 

jj 

.000 

.970 

3.   ; 

.157 

.001 

.951 

.132 

.001 

.953 

.116 

.000 

.960 

0 

i.  1 

.107 

.000 

.961 

.090 

.000 

.977 

.080 

j 

.000 

.963 

2.  ! 

.092 

.000 

.951 

.080 

.000 

.951 

.071 

> 

;ooo 

.955 

3.  \ 

.093 

.000 

.954 

.081 

.000 

.945 

i 

'   .071 

j 

.000 

.958 

1.2 

i.   ! 

.234 

.002 

.978 

.200 

.001 

.963 

.176 

i 

.001 

.964 

2-   1 

.213 

.001 

.961 

.183 

.001 

.958 

.163 

i 

.000 

.952 

3. 

.217 

.  001  • 

.954 

.186 

.001 

.961 

j   .165 

j 

.000 

.945 

82 

i.    ! 

.943 

.036 

.966 

.799 

.017 

.959 

.703 

i 

.011 

.961 

2.    i 

.855 

.023   ' 

.942 

.732 

.012 

.948 

.647 

\ 

i 

.007  I 

.945 

3. 

.871 

.026 

.936 

.743 

.012 

• 
.954 

.657 

i 

.008 

.936 

11 

i.    i 

4.54 

.822 

.954   ! 
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AH  =  A 
D020 
SH=S 
AVH  = 
VAR1 
VAR2 
VAR3 
CONT 
PR  =  M 
PC  =  N 
PO=N 
SH=S 
VAR1 
VAR2 
VAR3 
AVH  = 
WRIT 
FORM 
ITY«  ) 
WRIT 
FORM 
WRIT 
FORM 


FT  A 

-K  + 

-1 

XVA 

ET) 

INU 

L(K 

RET 

L(  I 

MOD 

=  L0 

INU 

)=A 

BIN 

CUN 

B+  1 

=  0 

T  =  r 

INU 
LRF 

=  0o 

=  0o 

=  0o 


L (K ) )305, 300,305 
2 

L(PRCBL,  Al,  Rl) 

2C1, 302,302 

E 

)=BET 

AL(  K) 

)=-ALnG(la-PLP)/T 

L ( I )-AMDA)310,310,311 

UN+1 

c 

MPD(  I  )-AMDOL( I ) 

+  BB(  I  ) 

+L0UNT-2) 175,174,175 


D/REP 

RIN/RF 

H/  R  E  P 

=SCAM/ 

H/LREP 

J  =  1  ,  LR 

H+(C AM 

AVH+(  A 

=VARU 

=VAR?+ 

=VAR3+ 

IMUE 

B/REP 

C  /  R  F  P 

F7REP 

H/(REP 

=  VAR1./ 

=VAR2/ 

=VAR3/ 

AVH/(R 

E(6,21 

AT  (  1  X  , 


D 

REP 

PP 

0A( J)-SCAM)**2 
HAT( J)-AH)**? 
(0(  J)-OD)  -<:2 
( PB( J)-BIN)**2 
(CC< J) -CM)** 2 


-U  ) 

(RFP-lo 

(REP-1 ) 

<REP-1, 

EP-lo) 

) 


3X,»SAMPLE    SIZE    LAMDA    CUTOFF    TIME    PRO 


E(6,22)    LRFP,N, AMOA,T, PROR 

AT(//,2I8,lX,3F10o5,//) 

E(6,23) 

AT(14Xt«n0Do     PROC  ,7X,  'BIN,     METH • , 6X , f CH I -  SO,     METH 


D,RIN,CH,VAR1, VAR2 ,VAR3,Pn,PB,PC 
VEPAGE' ,3( 5X,F 10 o5 ),//,! X, 'VARIANCE 


•,3(5X, 


WRITF( ^,26)0 

FORMATdX,  •  A 
)  t //, 
11X,  'COVERAGE  PROR1  ,3(5X,F10->5)) 

WRITF(6 ,24 ) AH, AVU 

WRITE(6,24)  SCAM,SH 

FORMATt IX, • AVGo  LAMDA  HAT= • , Fl 2o 6, • VAR I ANC E  LAMOA  HAT= 
,6) 

CONTINUE 

FNO 
ITERATES  TOWARD  LAMDA  FOR  PROCEDURE  1 


FUNCTION  XVAL(  ALPHA, Al  ,BU 

DATA  EPS/o5E-07/, N/1000/ 

A  =  A1 

R-Bl 

AL=ALPHA 

XLAST=-ol 
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M  =  0 

X  =  .l 

DX=„  1 
1.1  CALL  BOTP(  X,A,B,P,  P,  IFR) 

TEST=P-AL 

rF( TEST) 13, 5,13 
13  IF(APS(  TFST|-o^-O^I  12,12,6 
12  IF (A BS<  (TEST-XL  AST ) /XL  A  ST) -EPS ) 5, 5, 6 

5  XVAL=X 
RFTUPM 

6  IF(TESr*XLAST)7,8,3 

7  DX=-o5*DX 

8  X=X+DX 
XLAST=TFST 
M  =  M+1 
IF(M-N)11,9,9 

9  WRITE (6,1 5) 

15  FORMAT* IX, 'NO  CONVERGENCE1) 
XVAL=-U 
RETURN 
END 
CVAL  ITERATES  TOWARD  LAMDA  FOR  PROCEDURE  3 


11 


13 

12 

5 

6 
7 
8 
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FUNCTION    CVAL( ALPHA, G) 

DATA    EPS/, 5F-07/, N/IOOO/ 

AL=ALPHA 

XLAST=-ol 

M  =  0 

X  =  ,l 

DX  =  „1 

CALL    COTR( X,G,P,D, IER) 

TFST=P-AL 

IF(TEST) 13,5,13 

IF(ABS(TEST  )-o5F-04)     12,12,6 

IF(ABS(  (TEST- XL A  ST) /XLAST)-FPS)5, 

CVAL^X 

RFTURN 

IF(TEST*XLAST)7,8, 8 

DX=-o5aDX 

X=X+DX 

XLAST=TEST 

M=M+1 

IF(M-N)11,9,9 

WRITE (6,15) 

FORMATdX,  'NO    CONVERGENCE') 

CVAL=-U 

RETURN 

END 


5,6 
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